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ABSTRACT 

calihraH Jnn^ 31 challenge ® of th ® MSX science instruments require careful characterization and 
calibration of these sensors for analysis of surveillance experiment data. Procedures for reduction of 

^ (B f°> wm ■» praented is „° f f 

Pxnmnwt radiometnc Cand photometric) data and calculation of a precise MSX ephemeris. 
fn TaranJ i £ ® upport reduction ’ and the8e taken from ground-test data similar 

demin«r t ^ 1CS i° w MS3 S 8ens0rs and fr0m the IRAS satellite RSO detections. Examples to 
demonstrate the calculation of a precise ephemeris will be provided from satellites in similar orbits 
which are equipped with S-band transponders. 

1.0 INTRODUCTION 

The Midcourse Space Experiment (MSX) is scheduled for launch in 1994 and is a space Right 
program designed, m part to demonstrate surveillance of the space and Resident Space Object (RSO) 

tafraS™ n8o7(SPIRIT^rn t ’’ 11,6 ‘ e ' hni 7 challen e K! include cryogenic technology for cooling the 
infrared sensor (SPIRIT III), low noise high performance focal planes, high off-axis stray-light rejection 

be 898 km°Iltit T™ proces ® mg aad data compression, and contamination control. The orbit is specified to 
be 898 km altitude circular at nearly a Sun-synchronous inclination of 99.16 degrees. The lifetime of SPIRIT 

nerie/nf fin d ^ J lm0nths >. and the vis ! ble < SBV ) and UV (UVISI) sensors have a planned operation 
L mon ?• Prunar y saence data will be stored on board using tape recorders and downlinked via 

Srtl C rr Um . C l tl0n /V and com f ressed science data wiU be downlinked via 1 Mbit/s communications. 

w be A 1 l nd transponders on board > which allow the S-Band Ground-Link Stations (SGLS) 
network to provide tracking data for precise ephemeris determination. 

„ _ is tbe Principal space surveillance sensor and uses a 15 cm aperture off-axis re-imaging 

le ™^' a , th r m<KelCCt 7^ C CD f<*al P^e, a signal processo“nd sup^SS 

electronics. The SBV focal plane contains four three side abuttable frame transfer CCDs with 420x420 27um 
pixels each. The design characteristics are given in Table 1. ’ ^ Un 


Spectral Range 


Spatial resolution 


Field of View 


Aperture, f/no. 


FPS size (four CCDs) 


Frame times 


Quantum efficiency 


Table 1:SBV Characteristics 


0.3000-0.9000 um 


12.1 arcsec (60 prad) 


1.4° by 6.6* 


15 cm, f/3 


420 by 1680 pixels 


0.4,0.5,0.625,1.0,1.6,3.125 sec. 


28% 



The SPatial Infrared Imaging Telescope (SPIRIT) 3 sensor is the primary instrument on MSX 
covering the spectrum from the midwave infrared (MWIR) to the very-longwave infrared (VLWIR). SPIRIT 
III consists of an off axis re-imaging telescope with a 35-cm diameter unobscured aperture, a six-channel 
Fourier transform spectrometer, a five-band scanning radiometer, and a cryogenic dewar/heat exchanger. 
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The sensitivity of the spectrometer limits its use for observing RSOs. The radiometer has five Si:As focal 
rianTay^fi 8x192 pixels each, operating between 11 and 12 degrees Kelvin. It collects data m six color 
bands with a spatial resolution of 90 pradians. The scan mirror can remain fixed or can operate at a constant 
0 46 degree/sec scan rate with programmable scan fields of 1x0.75, 1x1.15 and 1x3.0 degrees. e ra lome r 

and Data Processor (OSDP) for clutter rejection and data compression. 


Table 2: SPIRIT III Radiometer Passbands 


Radiometer Bands I Band A I Band B1 I * ana v, I Dana - I — - 

I I RnnH B2 I I I 


6.0-10.9 



Band B1 
Band B2 


4.22-4.36 

4.24-4.45 


2 


10 


Band C 


Band D 


BandE 


11.1-13.2 I 13.5-16.0 I 18.1-26.0 



Passband (|im) 


Active Columns 


Sensitivity (NEFD) 

(10- 18 W/cm 2 ) 


The Ultraviolet/Visible Imaging and Spectrographic Imaging (UVISI) sensor system consists of five 
. nm-anliic imatrprs (SPIMS) and four imagers. Together the SPIMS cover a spectral range fro 
SSSEoSSSToS SSlsZm). Ae imagers include wide fold-uf-view (WFOV) and narrow- 
field-of-view (NFOV) sensors in the visible and ultraviolet. Surveillance investigations will co ^ entrate o ^ 
the NFOV imagers The commendable filter wheel in each imager houses three bandpass filters and a 
neutS densXmter, in addition to an "open” and "closed" position. The UVISI imager characteristics are 

given in Table 3. 


Table 3: UVISI Imager Characteristics 


Instrument 


FOV (deg) 


Resolution (p rad) 


Passbands (nm) 


open 


closed 


ND filter 


WB1 filter 


WB2 filter 


WB3 filter 


UV NFOV 


1.28x1.59 


90 



(xlO- 3 ) 


200-230 


230-260 


260-300(polarization) 


This report will focus on the reduction required for accurate analysis of the MSX space surveillance 
data. Three primary areas are identified: metric calibration, photometric and radiometric calibration, and 
calculation of the MSX precision ephemeris. 

The SBV and SPIRIT III should provide metric measurements in the FK5 reference frame accurate 
to 15-20^micr^radi£ms (3-4 arcseconds), and calibration will involve reference RSOs with well-known orbits 
touch as^geos and EGP) and calibration of the senaor boresite. MSX fiducial [reference frame, and censor 
alignments. The UVISI sensors are not expected to provide high quality metric data. 

The MSX sensors are designed to provide high quality radiometric data. The band- to band ratios for 
Vhmted SPIRIT III iiTadiance measurements are expected to be accurate to 5%. Radiometric call ra i 
S fSSS. San" 1 sMlarXence sources and multi-spectra! comparisons from the different sensors. 
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The reduction of SGLS tracking data should provide an MSX orbit which is accurate to 15 meters. 
Obtaining this accuracy is difficult, however, and involves calibration of the SGLS tracking data, modeling of 
the MSX attitude and cryogen flow rate effects, and modeling other non-gravitational effects. 

2.0 METRIC CALIBRATION 

The SBV and SPIRIT III sensors should be able to provide good quality metric data. It is expected 
that the SBV metric data will be accurate to 4 arcseconds or better and that the SPIRIT III will be only 
slightly worse. This section will first describe the reduction that is required to produce the most accurate 
metric observations, and then techniques of calibrating the data will be presented involving the use of 
precise orbits of calibration RSOs. Test data similar to the SBV will be shown as examples of the reduction 
procedures and the means by which these data are calibrated. 

2.1 Reduction of SBV Metric Data 

The SBV involves highly distorted optics (which are not defraction limited) due to a design which 
attempts to maximize the rejection of stray light from the focal plane. The size of each pixel is approximately 
13 arcseconds (60 micro-radians), and it is reasonable to try to sub-divide a pixel by a factor of 3-4 (or 4 
arcseconds) using light which spills into neighboring pixels. SBV is a self-calibrating instrument, since it will 
be able to detect stars down to 15 th visual magnitude. Therefore, a precise inertial location of the boresite 
and attitude map can be determined for the SBV at any instance by matching star detections to a star 
catalog and fitting attitude model parameters to the star map. Then, once an attitude model is updated for 
the SBV, the model can be inverted to focal plane locations for an RSO detection to produce right ascension 
and declination measurements. A metric error budget for the SBV involves accuracy of the reference star 
positions and number of stars observed, centroid error from the observed stars in SBV focal plane 
coordinates, (streak) endpoint error from the observed RSO detections in SBV focal plane coordinates, and 
model error from the SBV distortion map. 

Any catalog of reference stars can be used for SBV metric reduction, and we have chosen several. 
Each reference star is to be anchored to the FK5 inertial reference frame, and therefore stars which do not 
come in that frame must be carefully transformed. This is adequately detailed in (Smith et. al.). The catalogs 
which are currently implemented in the reduction software for the SBV are the SAO catalog, the 
Astrographic all-sky star catalog, the Guide Star catalog, and the Landolt Special Area Fields catalog (see 
Landolt). The Landolt catalog positions are not intended for metric use, but the color information is useful 
for photometric calibration of the SBV. Stars in the reference catalogs typically have position accuracy of 0.5 
- 2 arcseconds, and the use of dense catalogs often allows more than 100 observed stars from the SBV to be 
matched. Both annual and diurnal aberration are applied to transform mean to apparent place for each 
reference star, because diurnal aberration can produce an effect on the order of 5 arcseconds due to the 
velocity of the MSX about the Earth. The parametric model for the SBV attitude and distortion involves 37 
coefficients. 

The SBV will have a signal processor on board the MSX, which will detect linear streaks moving 
from frame-to-frame across the focal plane and stationary light sources. If the SBV is commanded to observe 
in a Sidereal mode, then the stars in the background will be stationary and streaks may represent RSO 
targets. If the SBV is commanded to track a particular RSO, so that the RSO is stationary in the SBV focal 
plane, then stars will be seen as streaking across the focal plane. Because the SBV optics are not detraction 
limited, a stationary point source will spill into neighboring pixels and a streak will actually appear as a 
swath of pixels between 3 and 5 pixels wide. This works to our advantage, however, since centroiding will 
allow the reduction software to sub-divide a pixel by a factor of three or four. Based on ground calibration of 
the SBV sensor, it is estimated that centroid error is less than 0.2 pixels (or 2.6 arcseconds). 

The direct map between inertial coordinates and pixel coordinates on a particular CCD in the SBV 
focal plane is accomplished by a series of coordinate frame rotations about angles which characterize the 
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SBV attitude. If v denotes an inertial unit vector and w a unit vector of direction cosines on the focal plane, 
the direct map is described mathematically as 

w = Ri(-p)R2(^)R3('I'+90)R2(90-5o)R3(ao)v, 

where Rk denotes the rotation about the axis k=l,2,3 by the indicated angle, ao denotes the SBV boresite 
right ascension, 5q the boresite declination, 4' the SBV roll angle, and p and £ denote off-axis angles which 
are designed to move the boresite vector from the center of the SBV focal plane to the center of a particular 
CCD. Once w is calculated, direction cosines are input to the distortion map so as to convert (linear) 
direction cosines to CCD pixel coordinates. 

The distortion map is expressed by a low-order polynomial in two direction cosine variables, U r and 
Uy, involving 32 coefficients. The map can be described mathematically by 

X c = 2 •muf m) lfl, <m) 

m=0 z y 

and 

= i 5 

m=0 

» 

where mathematical expressions for p(m) and q(m) are 
p(m) = [ m / 4 ] 

and 

q(m) = m - 4 p(m) 

with [ k ] being the integer portion of the number k. Initial calibration of the distortion map coefficients has 
been performed using the SBV optics on the ground, and it is found that this map is accurate to better than 
0.15 pixels (or 2 arcseconds). The map described above can be inverted to transform pixel coordinates to 
inertial space, but the distortion map must be inverted using iterative techniques due to its non-linear 
nature. 


2.2 Reduction of SPIRIT HI Metric Data 

In general, the SPIRIT III sensor will not be able to see enough stars in order to self-calibrate the 
pointing and attitude. Therefore, metric calibration will rely on alignment matrices between the SPIRIT III 
focal plane and an MSX fiducial frame. Alternatively alignment with the SBV focal plane can be used to 
determine the attitude of the SPIRIT III. The MSX is a rigid-body spacecraft, and there are no gimbaled 
mirrors. Therefore, pointing a sensor implies that the entire spacecraft must be pointed. The SPIRIT III and 
SBV sensors are co-aligned, and there exists a star camera on board the MSX which will be used to anchor 
the MSX fiducial frame. Whenever recognized stars are observed by the SPIRIT III, UVISI, and SBV 
sensors, then these data will be input to an algorithm to keep track of the alignment between sensors. This 
information will be steadily maintained throughout the lifetime of the instruments, and knowledge of the 
pointing for one sensor can be used in conjunction with the alignment information to determine the pointing 

for the other. 

Distortion in the SPIRIT III optics will be modeled in a similar fashion to the SBV, and it is expected 
for the coefficients to change little from ground-calibration values since the temperature of the sensor will be 
kept nearly constant. Star data will be used periodically to check and update the SPIRIT III distortion 

model. 

It remains to indicate how precise focal plane measurements of the RSO data are reduced from the 
SPIRIT III data. This task is performed by an On-Board Signal and Data Processor (or OSDP), which is to be 
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flown on board the MSX and can also be run on the ground with SPIRIT III raw data. The OSDP is a 
* re and _h ar ^ W y e system developed by the Hughes Aerospace Corporation, and involves two separate 
procedures. The first is known as the Time Dependent Processor, (TDP) which lines up the column data 
correctly from the SPIRIT III mirror scans (if necessary) and allows hot pixels to be grouped together The 
second is known as the Object Dependent Processor (ODP), which detects a group of hot pixels and identifies 
star and streak data. These object sighting messages are prepared and attitude models are inverted for 
conversion into an inertial coordinate frame similar to the SBV detections. 

2.3 RSO Calibration Orbits 

Precision orbits of calibration RSOs are to be used to check the accuracy of the SBV and SPIRIT III 
metric data and to determine biases for the data. This section will identify RSOs that can be used for this 
procedure and will describe the independent data sets that are used to calculate the calibration orbits. 

Good ^calibration RSOs are those which have stable orbits which are not difficult to model and which are 
routinely tracked by a variety of sites so that a dense sampling of tracking data is available. The Lageos 
calibration spheres are excellent candidates, because they are equipped with laser cube corner reflectors and 
serve as calibration RSOs for both radar and optical sites around the Globe. The Lageos orbits are stable and 
are routinely known to within 10 cm. They can also be routinely tracked by MSX sensors since viewing 
angles can be easily found which do not require tracking too close to the Sun or the Earth’s limb. Other 
RSOs which make good calibration orbits are EGP, Etalon, ERS-1, and GPS. The GPS satellites make good 

calibration orbits because accurate samples of state vectors accurate to better than 5 m are available at all 
times. 


The objective for calibrating the SBV and SPIRIT III data is to obtain an independent reference orbit 
or an o served calibration RSO which is accurate to better than 1 arcsecond in sensor right ascension and 
declination angles. The reference orbit is calculated using a special perturbation orbit determination 
program which can fit many types of data to the equations of motion using a detailed force model. Then this 
orbit is compared to the observed data and statistics are calculated from the pass. Biases are included in a 
historical database and incorporated into the calibration models for the metric data reduction. 

Tk e Pj*® cision orbit determination software which is used for this procedure is known as DYNAMO 
and it has a history dating back to the 1960s for use in calibrating sensor data and providing precision 
reference orbits. The capabilities and qualities of DYNAMO are highlighted in Figure 1 


Force Models 

• GEMT2 (42x42) Gravity 

• Mood & Sun Third Body 

• Atmosphere Drag (MS IS) 

• Cryogen Venting 

• Solar Radiation Pressure 

• Body Tides 

• Ocean Tides 

• Albedo Pressure 

• General Relativity 

• Thrusting 

| Reference Frames 

• Satellite Laser & 
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• J2000 & 1950.0 Inerti: 
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DYNAMO 
Accuracy lm 


[Tracking Date Types 
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] 


Weighted Least Squares 


Solve For 
Parameters 


• Radar 

• SGLS 

• Laser 

• Telescope 
•GPS 

• VLBI 


Atmospheric Drag 
Cryogen 
Geopotential 
Solar Radiation 
Station Coordinates 
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Thrust 

Station Biases & Re fraction 
Output Options 

• Propagated 

Earth-Axed Cartesian (DMA) 
ECI Cartesian 
Kepler Elements 
Station Topocentric 

• Extensive Diagnostics 

per Observation 


Figure 1: DYNAMO, A Precision Orbit Determination Software Package 
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2.4 Example of Metric Reduction and Calibration 

In this section, an example of the metric reduction and calibration process will be granted. The raw 
data for this example are selected from a ground-test setup designed to be similar to the SBV_ A 6 mch 
telescope system was attached to the 31 inch optical mount at the Experimental Test Slt ^™) located 
near Socorro, New Mexico. A 420 x 420 CCD focal plane was used to record the data employing parameters 
for the integration time and number of frames which are similar to the SBV. In November of 1991, a 32 
minute pass^f the Lageos I satellite was collected with this system ^ 29 fte of data with 8-16 frames for 
each set The reduction and calibration steps for this set of data will be highlighted m this section, with the 
final objective being to characterize the accuracy of Lageos I metric observations. 

First, the raw data must be passed through the SBV Signal Processor to detect the stationary point 
sources and streaks. In addition to the streak metric information, a type of signature data must be collected 
which labels frame and intensity data for each pixel near the best-fit line to the streak. The signa ure 
information is used in post-processing to refine the end points and deduce a visual magnitude measurement 

for each streak detection. 

The next step is to match as many of the observed stars as possible to an on-line star catalog and to 
update attitude parameters for the sensor (boresite vector and focal plane scale) The results are shown in 
Figure 2 displaying a Bull’s-eye plot of the residuals (in the focal plane) of the stellar position data after the 
fit This’ chart shows that the root mean square (rms) of the residuals is close to 1 arcsecond in eac 
direction. Outliers can be traced to saturation, double-star systems, or stars which are on the edge of the 

field of view. 


|rms error 

X — 1.10 arcsoc 
Y — 1.19 arcsec 




1 PIXEL 


I I I I I I H 



Figure 2: Bulls-Eye Plot of Residuals (arcseconds) of Matched Star Positions Based on a least Squares Fit of 
the Focal Plane Attitude Parameters 

Once an accurate attitude map is available, the end points of the observed streak must be refined 
using the signature data and the attitude map must be inverted to transform the streak end point 
measurements* in focal plane coordinates to inertial right ascension and declination measurements. In order 
to characterize these measurements, a precision ephemens for Lageos I must be de ^ erm _ 
accomplished using DYNAMO and independent measurements of Lageos I from a time period spanning pi 
or minus three days about the epoch of the pass. The independent measurements are taken from laser radar 
and skin tracking from radars and optical sites in the Space Surveillance Network. This suite of observations 
allows an ephemeris for Lageos I to be determined to better than 10 cm, and this im P h ® s that the pre * c £ 
senTor Accuracy from a ground-based site is better than a fraction of an arcsecond. A comparison of the 
Lageos I observations to the precision ephemeris is shown in Figure 3. There is a strong bias m declination o 
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2.2 arcseconds, and this could be due to a time discrepancy since the inclination of Lageos I (109 degrees) 
causes along track orbit error to be manifested as a declination error. The standard deviations in right 
ascension and declination are larger than the desired 4 arcseconds, and the principal sources of error are not 
easily identified. Possible error sources include a (slight) shift of the mount from frame to frame and noise in 
the data. 



RIGHT ASCENSION 


J? ulls Eye Plot Residuals (arcseconds) of Lageos I Metric Observations as Compared to a 
DYNAMO Reference Orbit for a 32 Minute Pass. 

3.0 RADIOMETRIC CALIBRATION 


3.1 Introduction 


The SBV, SPIRIT III, and UVISI sensors should be able to produce high quality photometric and 
radiometric data. Surveillance data will be used with high fidelity models of the reflected and self-emitted 
radiation from satellites, to develop methods for identifying RSOs, monitoring their status, and determining 
some of their physical properties. These high fidelity models account for such variables as: solar phase angle, 
sensor to RSO aspect angle, material properties, and temperature. The discussion here is limited to the 
methods of data reduction and calibration to be used with each of the MSX sensors, and an assessment of the 
expected accuracy. 


Each of the MSX sensors is unique, and the data reduction and calibration will be different for each 
one. The MSX program intends to provide data that is certified to be calibrated to within specified limits To 
this end the Data Certification And Technology Transfer (DCATT) Principal Investigator Team is devoted to 
establishing calibration procedures and standard data reduction software. Even though fundamental 
differences between the MSX instruments require different calibration and data reduction methods, they do 
share some common elements. For example, there are three steps in the calibration process Each 
instrument will have extensive preflight (or bench) calibration. It is hoped that this preflight calibration will 
be valid for the on-orbit data. Second, there will be a series of on-orbit observations taken for purposes of 
calibration. These include internal sources, used in each data set, to obtain corrections for each data set. 
Finally, some data sets will contain observations of objects with known luminosity, which will provide 
additional calibration information. These three sources of information will provide the calibration 
information necessary for photometric and radiometric analysis of MSX data. 
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3.2 SBV Radiometric Reduction and Calibration 

The conversion of raw digital numbers to engineering units for the SBV is done for eaC ^ p1 ^^ 
dark current (a function of temperature) is measured with the SBV cover closed and is subtracted from t 
measurement. The responsivity conversion of digital numbers to Watts/cm in band is done using groun - 

based calibration data. 

The on-orbit calibration phase will be done using calibrated reference stars primarily the Landolt 
fields (Landolt). The Landolt star fields have been established as astronomical photometric calibration 
standards Since the SBV has a nonstandard and very broad spectral response the relation between the S 
magnitude and the magnitude of a star must be known. From analysis of the SBV spectral response the 
”Enw“ee„ the visual magnitude of a .tar, M., with color, B-V, and the SBV magmtude. given 

by an empirical relation (Beavers) 

M sbv = M v + Ri(B-V) + R2(B-V)2 + R 3 (B-V)3. 

We expect this relation be have an accuracy better than 0.02 magnitude. 

During the collection of RSO science data, many observed stars will havea visual magnitude and a 
color though not of the accuracy and reliability as the Landolt calibration fields. The process of reducing the 
SBV d^U start with a computed value, M SBV - (as described above) for each detected star .found in th 
star catalogue, and this is compared with the observed magnitude, M SBV ° ba , derived from the calibration, 
mean difference between these two quantities becomes the zero-point correction for the frame set, 


AMsbV = 21 ( M SBV obs - ^SBV^ ^ n 

and is subtracted from the observed RSO magnitude. In this way the SBV makes self- calibrating 
photometric measurements. 

In processing SBV photometric data, a simple model is used for prediction of Msb^ 0 . It is based on 
assuming that the satellite reflects the solar spectrum. We have adopted (Beavers) 


Msbv" 80 = -26.8 ‘ lo Sio < P A /r2) 

where p is the reflectivity, A is the effective area, F(<t>) is the phase function, and r is the range to the target. 

A, a;, e*mnple of tin. pruces., a ^ pa» of i SBV Ito ££££> 

safTcxampl" datataMsseY^the Section 2.4 on metric calibration, andphotomctric measurOTents^ccre 

residuals in mv gm reflectivity A=0.5 m 2 for the effective area of the sphere, and F(<t>)=2((7i 

sphere. The results illustrate modeling the photometric properties of 

Lageos sample to better than 0.4 magnitudes. 

3.3 SPIRIT HI Radiometric Reduction and Calibration 

rm. „i:u ra t; n n nha<?e for SPIRIT III will be done in two different ways. First, a ground 
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a 8 .ub^ffll1“ P d “;T " f “" da '”“W "ay on the 

issue in the calibration oflRAS (Beichman SC£^ Wl T mant .\, Thls !»*«• to be an important 

Se , nS ° r ’ but another calibration method will be used^e^ StabIe reference to monitor a 
that will be deployed from the MSX during the SPIRTT TTT 1ST' five emisslve reference spheres 

with Martin Black, are designed to haTaweU define! 8pheres > ^ oated 

will sample the full dynamic range of the SPIRIT ID TT^ aSl i to Provide orbital geometry that 
reference sphere temperature with sufficient accuracy todet & ^ P ^ 0Vlde knowl fdge of the emissive 
band-to-band ratio to 5%. y 10 determme the absolute irradiance to 15% and the 


Table 4: 


OBJECT 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 
8820 


Comparison of an Observed Pass oQ «sible Data for Lageos with a Photomeirio Model 

SENSOR MAGNITUDES 


YR DAY 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 
91 312 


HR MN 
2 6 


2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


6 

6 

7 

8 
8 
9 
9 

15 

16 
18 
18 

19 

20 
21 
22 • 

23 ■ 

24 - 

25 - 

26 - 

27 - 

28 - 

37 - 

38 - 


PHASE 
+ 26.1 
+ 25.6 
+ 25.3 
+ 24.4 

- 22.7 

- 22.7 

- 21.6 

- 20.9 

- 17.9 

- 18.7 

- 20.7 

- 21.7 

- 22.9 

- 25.0 

• 27.0 

29.2 

• 31.6 

34.2 

36.3 

38.6 

41.2 

43.6 

64.2 

65.2 


RANGE 

7581.1 

7548.4 

7534.7 

7480.7 

7382.1 
7382.1 

7312.9 
7270.0 

6880.7 

6851.9 

6814.4 

6804.7 


6798.8 
6798 
6806 
6823 
6850 
6891.0 
6930.2 

6981.5 

7045.6 

7114.7 
8068.4 

8131.8 


OBS 
12.10 
11.98 
11.93 
12.10 
11.72 
11.58 
12.02 
12.13 
11.71 
11.92 
11.90 
11.88 
11.76 
11.55 

11.54 
12.06 
12.13 
12.53 
12.80 
12.67 
12.81 
12.37 

13.55 
13.24 


PRED 
12.21 
12.20 
12.19 
12.17 
12.13 
12.13 
12.10 
12.08 
11.95 
11.94 

11.94 

11.95 
11.95 
11.97 
11.99 
12.01 
12.04 
12.08 
12.11 
12.15 
12.20 
12.26 
12.86 
12.90 


RES 
- 0.11 
- 0.22 
- 0.26 
- 0.07 
- 0.41 
- 0.55 
- 0.08 
0.05 
- 0.24 
- 0.02 
- 0.04 
- 0.07 
- 0.19 
- 0.42 
- 0.45 
0.05 
0.09 
0.45 
0.69 
0.52 
0.61 
0.11 
0.69 
0.34 


quick check of ttaXereed RSoTadiXTlry^eXl 1 ' 81 ' 0 " a nSe ? d f ° r 016 automat « d Processing and a 
developed for processing the six SPIRIT HI wavebanX and fliTf* 8 Xao process - Tables have been 
augmented for other wavebands as necessary The " IRAS wav ebands. These can be 

each band in Watts/m 2 in band to any oSr h LTl.IT conversion of the input observation for 

(W/m 2 /Hz) and Naj’s (W/m 2 /micron) The baste rtSLtnfa L XT"” 7A* an,lard units such as Jansky’s 
alter + blocker * detector response for *?■“>* * WM, the 

absolute temperature, then we define ,be ob*c, £j£S radian't T “ «“ 

r 00 


f b (T) = J* a)3 X (T) dx 01) 


where the Planck flux density is 


3 X (T) = 


3.74185 x 10 8 


© 


H^-i) 

flux density unit little Naj. A generally used alternate 

microns/sec, we can convert Nafs to JanskX with EiT lhat f«=2.99792458xl0“ 

the observed flux (in-band) irrad ance as J s ' da nsky's*2.99792458xl ( H 2 A 2 We can now compute 
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f ’ m - m im£T) iX 

where A is the area of the object, e the emissivity, and r the range to the object. 

BuXj F ^.r P ssr t hi 3 bi rro n dt"^ the rati ° 

_ Fq(D = F*(D 

for bandstand y[is^ moiwtliflc function of T. ^ 

SPIRIT III sensor. emissivity. With the measurement of n m band Runes. 



1 ®fnpor»iui» v»v 

Figure 4: Ratio of Flux Densities for SPIRIT III Wave Bands 
^ U3f Given a temperature, the observed in-band flux can be expressed in Naj’ 8 at a reference wavelength, 


V 


as 


0 - °V F * ( T) ’ ^ h band This can be converted to Jansky’s, as 

t is customary (though not necessary) to select \ m 
[escribed above, with 


•=N?( 


i 


2.99792458 x 10 


, 12 ) 


, - f Vio traS satellite. The InfraRed Astronomy 

As an example of this process, we show some data ^ ^ IRAS sensit ivity is very similar to that 
’’atpllite (IRAS) was operational for about 10 mon cronmptrv resulted in measurements at a phase 

expected from SPIRIT III, although the on medium to high altitude satellites were made 

“ f ei"t y ua s r e o Lr a 

<» -* ° nhe i2 'tft 0 a *e 

ssasa: sar a* ^ =— — to be hi8h 11 is reasonab 


72 


methodology to be able to provide temperature data from SPIRIT 
and ±2% in emissivity and absorbtivity. 


Ill measurements to ±5*K in temperature 


3.4 UVISI Radiometric Reduction and Calibration 


issues SESE2L1! H ?T ““rated »•- - — . There are two 

^ apd ihe ™ si — 

Table 5: Comparison of an Observed Pass of IRAS Infrared Data for LCS-1 with a Radiometric Model 

OBJECT YR DAY HR MN END RANGE 5? (W/SQ mBtL) 

B s I a 1 1 El “1 fl || ifs ^ & 

3 TEMPERATURE MEASUREMENTS: 479.51 561/76 9 5 q 9 7 J- 11E " 13 °' 14E - 13 -°- 27 E-14 

erence between these two quantities becomes zero-point correction for the frame, 

AMuvisi = ^ ( MuvjgjObs . Muscat )/ n 

phetomehicmeMui^Mnto)^imaar'totI^^^ ad11 *^ Ude * Ws Way ** ™ SI —• «* calibrating 

In processing UVISI photometric data, a simple model is used for prediction of M Rso Tf ■ . 0 , 

on assuming that the satellite reflects the selar spectrum. We have adopteS?Beave“ ) ^ ^ 

M UVISI RS0 = - 26 -8 * log 10 ( pA F(<J>) /r 2 ) + AM , 

SSSrS “ d ' is tha *■**• 

4.0 THE MSX EPHEMERIS 


4.1 


Description of the Problem and Method of Attack 


the Z th r PhemeriS " P ° Sm ° " accuracy of 

how it can be achieved. d * re< ' u,r,!d epherncns accuracy for the MSX and 


SPIRIT ni C a e „ P d h SRVTt a “ UraC ^ r ^ r ? m “ U are dep<!nde ' lt on the required metric data quality for the 

asas^ssai 
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that with 4-arcsecond sensor data quality and a requirement that the ephemeris error ba at leant 3 .mailer 
implies that the ephemeris accuracy for MSX must be better than me rs. 

Thera ware three method, considered 

GPS receiver, 2) ground (or skin) tracking from the Spa receiver on board the MSX was the 

tracking from the Air For* S-band ; Ground^ ^Af^^quirementa. Radar (skin) 

most attractive choice, but it was elimina Millstone radar in Massachusetts and from the 

tracking of the required These provide accurate enough 

ALTAIR radar in the Marshall Islands (in the near ^uahmtu ramlicy ^ ^ ^ 

measurements, but they are heavily taske a rea >' an network is used by the Air Force for satellite 

MSX would demand too much of their resources. Th : SGLS network Msx ^ 

communications, and the S-band tracking a aao SGLS network Tracking data from this 

have a coherent S-band transponder and wfl h > Se MSX. Thf SGLS network 

— , and elevation for 

satellites with a coherent transponder. 

To determine if the SGLS tracking data could slu^etwere™! in 1^1989^ 

MSX, a number of evaluations of the data h^eb^n mad . TCi similar Qrbital paramete rs to 

through 1990. These studies involved SGLS data from two sateu paramete rs for MSX. The two 

MSX. Table 6 compares these with one c0 ^ n de 7* e f ^„ th ® r ^ not well-determined, which implied that 
test satellites were lower, however, an eir : a ^ , g ected for MSX. The remaining dynamical 

atmospheric drag was more difficu fact that MSX will have cryogen 

jrJSSS Sir m sensei r. Phenomenon will constitute a 

significant perturbation of the MSX orbit and is absent in the test satellites. 


O; Utrir 

MSX 

19911 

20497 

ALTITUDE (km) 1 

ECCENTRICITY 
INCLINATION (deg) 

DRAG 
MANEUVER 
SATELLITE ASPECT 

888. 

0.001 

99. 

YES 

NO 

KNOWN 

490. 

0.00133 
47.7 
HIGH 
OFTEN 
NOT KNOWN 

460. 

0.00161 

43.1 

HIGH 

NO 

NOT KNOWN 


The objectives ofthe of SGLsS i^ metersVr^lge^ 

mdlidegrees for the^^— 

: v Tut^ - d what a “ 

processing has to be done in order to use the data to its potential. 

The late 1989 and 1990 studies took place during periods of major solar activity ^the^. 1 ^ 
altitude satellites, it was found that three 1 d^ntervals were also estimated. The orbit 

solving for the satellite radL measurements and also by comparing overlapping orbits^ 

“ar”me— L TUSSm the MUlstone H.il L-Band radar (with a range accuracy of 1 meter, and 
from the ALTAIR UHF radar (with a range accuracy of 10 meters). 

The SGLS data had a nominal range bias correction applied (none jangles, and a t = here 

^apphed'and so we applied corrections based on the best 
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available global model (Bilitza, Klobuchar). The range rate had to be converted to a range difference 
measurement to be properly used. Additional biases in the range data were determined thereby improving 
the calibration and accuracy of the data. The angles were also calibrated. 

4 ,,-x W ! th an e , nhanced calibration and an average of 8 tracks per day, the evaluation showed that the two 
satellites have orbits computed to an accuracy of 15 meters. This is illustrated in Figures 6 and 7. Figure 6 
provides a sample of the data quality from one of the SGLS tracking stations (Guam) using an orbit 
computed with the radar data, and Figure 7 shows range residuals for the radar data using an orbit 
computed with the SGLS data. The plots display the mean and +/- 1 a error bars for each of the tracks. The 
means from Figure 7 are within 15 meters and are due to a combination of measurement error and orbit 

^ ese evaluations indicate that the SGLS network can provide the necessary tracking data for the 
MSX orbit computation. 


SAMPLE ACCURACY OF DATA (Guam) 
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Figure 6: Residuals of SGLS Tracks from an Orbit Determined with Radar Data 


SGLS ORBIT ACCURACY EVALUATED USING 
RADAR RANGE MEASUREMENTS 



Figure 7: Residuals of Radar Tracks from an Orbit Determined with SGLS Data 
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4.2 Modeling the MSX Orbit 

In order to produce a 15 meter orbit for the MSX on a routine basis, accurate modeling of the forces 
acting on the MSX orbit is required. The gravitational forces are well known and already incorporated m the 
DYNAMO software. The non-conservative forces which are more difficult to mrid indude radiation 
pressure, atmospheric drag, and the force induced by the venting of the gas created by the sublimation of a 
solid hydrogen block within the cryo-stat. Two important data types that must be obtained and input to the 
MSX force model are the satellite attitude and the venting flow rate of hydrogen gas from the cryogen 
cooling system. The attitude data is necessary for the drag, solar radiation pressure, and cryogen venting 
models. The cryogen flow rate data will be measured and provided by Utah State University (the 
manufacturers of SPIRIT III), and details are still being worked out. 

The cryogen gas venting model is the remaining dynamical complication for the MSX satellite. The 
cryo-system on the MSX contains a large mass of solid hydrogen, which is used to keep the instruments at 
the required temperature. As heat is added to or generated by the satellite, the hydrogen escaping through 
the venting system produces a low thrust on the satellite. Depending on the detailed geometry of the venting 
system, the effect of the resulting force can be large when integrated over the course of a day. Details on the 
exact nature of this force are yet unclear, but an exact or an empirical model (which will be parametrized 
and updated from the orbit fits) is critical to achieving a 15 meter orbit for the MSX on a routine basis. 

5.0 SUMMARY 

The MSX science instruments cover a wide range of the spectrum from the ultra-violet to the long- 
wave infrared and will be able to provide useful metric, photometric, and radiometric data for surveillance of 
the Resident Space Object background. Careful characterization and calibration of the sensor data is 
required for accurate analysis of space surveillance experiment data. Procedures for reduction and 
refinement of the metric and radiometric (or photometric) data have been presented and methods of 
calibration have been described. Examples from ground-test data similar in characteristics to the MSX 
sensors and from the IRAS RSO detections have been presented to support the reduction and procedures 
outlined in this report. In addition to the reduction of the MSX sensor data, it is crucial to calculate a precise 
MSX ephemeris. The ephemeris is calculated using SGLS tracking of data from one of two - an 
transponders on board the MSX and sophisticated models of the MSX orbit. Examples to demonstrate the 
techniques for this calculation have been provided from satellites in similar orbits equipped with S-ban 

transponders. 
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